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ABSTRACT 



Complex temperature anomalies observed near Arctic 
ice margins in the Chukchi Sea were found to be associated 
with the interaction of the warm coastal current with the 
ice cover and resident bottom water. These anomalies were 
characterized by large temperature gradients found in areas 
of very low density gradient. They were in the form of 
inversions, interleavings, and discrete parcels which varied 
greatly in short distances. The processes which produced 
mesostructure were directly linked to the presence of ice 
and were found to affect the entire water column. Structure 
was found in the vicinity of the ice margin but was observed 
to dissipate well inside the ice margin. A dynamic high was 
found which was related to the melting of ice and character- 
istically occurred in the melting zone of the ice margin. 

4 

This dynamic high was limited to a depth of 10 meters and 
is believed to be a factor in the deepening of mesostructure 
elements and other phenomena as they pass the ice margin 
with the current. 
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I. INTRODUCTION 



A. DESCRIPTION OF THE PROBLEM 

An investigation was conducted to examine the mechanisms 
associated with the phenomenon of temperature mesostructure 
found in the water columns near Arctic ice margins. The 
term mesostructure describes the relatively large anomalies 
which occur in the vertical temperature profiles; such a 
profile is shown in Figure 1. Corse (1974) defined an 
element of mesostructure as that portion of the profile 
from the depth where the temperature gradient became positive 
to the depth at which the temperature was again the value 
it had when the gradient became positive. These anomalies 
were considered as warm, vice cold, since they were found 
to be associated with the advection of warm water into the 
region. The magnitude of an element was found to be of 
the order of meters in depth and of tenths of degrees in 
temperature variation. The phenomenon of mesostructure 
was found to be associated with the presence of an ice 
margin and was present to some degree in all ice margin 
crossings investigated in this study. 

B. DATA BACKGROUND AND SOURCE 

Investigations into temperature structure near ice 
margins had their origin in reports of severe deterioration 
of sonar propagation, which was believed to be the result 
of complex sound-speed profiles near the ice margins. As 
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Figure 1. Definition of Mesos tructure Elements. 



a result of this interest, research cruises, under the 
general project title MIZPAC, were directed by the Arctic 
Submarine Laboratory, Naval Undersea Center. These cruises, 
carried out in 1971, 1972, and 1974 on Coast Guard ice 
breakers, have resulted in a data base of 350 STD stations 
with another such cruise scheduled for 1975. The results 
of the 1971 and 1972 cruises were published (Paquette and 
Bourke , 1973, 1974) and those for 1974 and 1975 are in 
preparation by Paquette and Bourke. The data from the 1974 
cruise were analyzed because they were the most useful of 
the available data, due to the number of ice crossings and 
the close proximity of the stations. 

The 1974 data were collected using a Bissett-Berman 
Model 9006 STD. Two lowerings were required for each 
station; a deep cast, and a shallow cast with the conductivity 
sensor shunted with a precision resistor to enable recording 
of the low salinities encountered near the surface. Seven 
series of closely spaced stations were analyzed, each one 
representing a penetration of the ice margin and termed a 
"crossing." Figure 2 indicates the location of the crossings, 
carried out during the period 18 July to 26 July 1974, and 
several positions of the ice margin for this period. The 
positions of the ice margin were derived from cruise data 
and NOAA-4 VHRR satellite photographs interpreted by the 
NOAA-NESS Environmental Products Group. 
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Figure 2. Location of Crossings and Ice Margin. 



C. DESCRIPTION OF THE MIZPAC 74 AREA 



The Chukchi Sea is a marginal sea with very little 
bottom relief in the area of interest and with an average 
depth of about 45 m. 

The ice coverage and thus the ice margin varied 
considerably during the months of June through August 
1974 (Fig. 2). There is much year-to-year variation. In 
1974 the ice melted back from Bering Strait in June to 
73°N by the end of August, farther north than normal. 

The pattern of the melting process followed closely 
that of the previously known current system of the Chukchi 
Sea. The summer surface circulation of the Chukchi Sea is 
dominated by a general northward flow of warm water from 
Bering Strait. This flow was observed by Paquette and 
Bourke (1974) and was typified by a surface temperature of 
5° to 11°C during late June through September and probably 
involved the upper mixed layer of 10 to 20 m depth. Little 
is known about the movement of the lower layer. The current 
was described as staying to the right side of the Chukchi 
Sea and splitting at Cape Lisburne, one branch following 
the Alaskan coast and the other turning north-north -west 
toward Herald Shoal. Garrison et a. (1973) reported the 
speed, relative to an ice floe, of the current along the 
Alaskan coast near 72°N to vary between 10 and 50 cm sec \ 
with a velocity shear between 10 and 20 m which was usually 
associated with the pynocline. 
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D. OBJECT OF THE STUDY 



The thrust of this undertaking was to gain an under- 
standing of the mechanisms by which mesostructure is formed, 
utilizing new data and previous findings. It was particu- 
larly directed toward finding the relationship between the 
formation of mesostructure and the following processes: 
the influx of warm water and the melting of the ice cover, 
the sloping of the density surfaces, the dynamic heights 
near the ice margin, and local currents. Additionally, 
it was desired to determine if different processes were 
responsible for deep and shallow mesostructure. 
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II. INVESTIGATIONAL PROCEDURE 



The data presentations were predicated on the assumption 
that warm water flowed northward toward the ice and that the 
interaction of ice and water resulted in mesostructure . It 
was expected that the phenomena involved would be evident 
in vertical sections of the conventional oceanographic 
properties along the lines of the several ice-margin 
crossings. Temperature and sigma- t sections were therefore 
prepared. The latter, being closely related to salinity 
at locally observed water temperatures, eliminated the 
need for salinity sections. 

Later it became apparent that the question of direction 
of flow of the warm water could be investigated by means of 
the heat and ice melt contents of the water column. The 
local dynamic height gradient near the ice margin seemed 

< 

a possible driving force for downward mixing. These 
parameters also were computed for the ice-margin crossings. 
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III. DISCUSSION 



A. TEMPERATURE CHARACTERISTICS 

The temperature profiles, nested in sequence, for each 
ice margin crossing are shown in Figures 3 through 9. The 
station number and bottom temperature (°C) are indicated 
at the bottom of each trace, and the ice concentration in 
oktas (eighths) is displayed at the top with the symbol 
"X" indicating no data. The station numbers increased as 
the ice margin was penetrated from outside the ice; thus 
the crossings proceeded from left to right in all figures. 

The temperature profiles displayed several characteristics 
which were common to all the crossings in general. A layer 
of relatively warm water, between the surface and 20 m depth, 
appeared outside the ice in all cases. This layer was 
cooled and modified with progress toward and across the ice 
margin and was greatly diminished or absent well inside the 
ice margin. Surficial cooling caused the temperature 
maximum of this layer to descend several meters when ice 
was encountered forming the "nose" feature described by 
Paquette and Bourke (1973) . 

Corse (1974) described mesos tructure as being divided 
into three classes: the nose feature, and shallow and deep 

structure. He defined shallow structure as having sigma- t 
values less than 25.5 and deep structure as that with 
greater values. He placed the sigma-t value of 25.5 at 
about 20 m depth; however, as will be seen, this 
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Numbers at the top of trace are ice concentrations; 
at the bottom of the trace the bottom temperature. 
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Numbers at the top of trace are ice concentrations 
at the bottom of the trace the bottom temperature. 
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Numbers at the top of trace are ice concentrations; 
at the bottom of the trace the bottom temperature. 
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Numbers at the top of trace are ice concentrations; 
at the bottom of the trace the bottom temperature. 
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at the bottom of the trace the bottom 
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Numbers at the top of trace are ice concentrations 
at the bottom of the trace the bottom temperature. 



(SU3JL3N) HJLd30 





bO 

C 

•H 

CO 

10 

O 

U 

U 



co 

0 



•H 

m 

o 

u 

Ph 

0 

p 

3 

<P 

cd 

P 

0 

P. 

e 

0 

H 



co 

0 

P 

3 

bo 

•H 

P-i 



25 



co • 

P 0 
O U 
*H 3 
4-> 4-> 
cd cd 
P P 
4-> 0 

£ a, 

0 s 

O 0 
C 4-) 

o 

o e 

o 

0 4-> 

U 4-> 
•H O 
X> 

0» 

P (D 

cd 

4-> 

0 

L) 0 
cd u 
P cd 

4-> P 
4-> 

m 

O 0 
X 
P<4-> 
o 

4-> m * 

o 

0 

^ £ 

4-> O 
4-1 
4-> 4-1 

cd o 

to 

P 0 

0 rP 

X> 4-1 

6 

3 4-1 
^ Cd 



investigation found the 26.0 surface at about 20 m. This 
depth was used to make a distinction between shallow and 
deep mesostructure . 

Shallow structure was present in varying degrees of 
severity in all crossings. Crossing 3 exhibited the most 
severe shallow structure as compared to the lesser intensi- 
ties found in Crossings 4, 5, and 6 and to Crossings 1 and 
2, which displayed very little shallow structure. The 
presence of shallow structure correlated well with the 
presence of the nose feature, which was present in all 
crossings except 1 and 2. Thus, the nose feature appeared 
to have a causal relationship with shallow structure. 

Deep structure was found to be severe in Crossings 1 
and 2 and very slight or absent in the other crossings. 

The presence of deep structure appeared to be a separate 
phenomenon from that of shallow structure. 

The temperature profiles displayed no apparent 
correlation between the presence of mesostructure and the 
maximum temperature present. This differed from the 
belief of Corse (1974) in the relationship between high 
temperatures and the presence of deep mesostructure. 
However, the formation and presence of structure, both 
shallow and deep, did correlate with the presence of an 
ice margin, as did the formation and descent of the nose 
feature. The mesostructure eventually almost disappeared 
at some point within the ice margin. This may indicate 
that the structure had been dissipated as it moved 
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"downstream" with a northerly flowing current past the ice 
edge or that there was little or no water flow toward the 
ice . 

Vertical cross-sections of temperature for the seven 
crossings are presented in Figures 10 through 16. Again 
the relatively warm surface layer, with a maximum tempera- 
ture of 3 to 6°C, may be observed south of the ice at 
depths of less than 20 m. This layer cooled drastically and 
lost most of its warmth in the vicinity of the ice margin. 
This was obviously due to a utilization of this heat for 
the melting process which will be discussed in a later 
section. Referring to the cross-sections, the isotherms 
were very irregular and displayed large gradients, many 
inversions, interleavings, and lenses. These irregularities 
were more pronounced in the vicinity of the ice margin and 
dissipated when well inside it, as has already been stated 
in the description of the profiles. Individual elements of 
mesos tructure usually could not be traced for more than two 
consecutive stations, a distance of 1 to 5 km. The maximum 
distance over which a mesostructure feature was traced 
was about 11 km. This is in reasonable agreement with the 
results of Corse (1974) . He found typical correlation 
distances of 1 to 2 km. Structure was found to exist for 
some time after the disappearance of ice from an area 
since it was found as far as 28 km south of the ice margin 
where the ice had apparently been recently. 
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Figure 11. Temperature Cross Section, Crossing 
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Figure 12. Temperature Cross Section, Crossing 
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Figure 13. Temperature Cross Section, Crossing 
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Figure 14. Temperature Cross Section, Crossing 
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Figure 15. Temperature Cross Section, Crossing 
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Figure 16 . Temperature Cross Section, Crossing 



Also well correlated with the ice margin was the presence 
of cold bottom water which was absent to the south of the 
margin. This type of water appeared fairly precisely at the 
margin and thickened as one progressed into the ice. It 
had a temperature of -1.6 to -1.7°C and was usually either 
in the form of a wedge, as in Crossings 1, 2, 3, and 7, or 
that of a cold tongue at about 30 m depth as in Crossings 4, 
5, and 6. This correlation of the boundary of the cold 
bottom water with the ice margin is indicative of the trans- 
fer of heat downward in the water column near the margin. 

The dissipation of the surface heat across the margins, 
the melting of the ice, and the absence of the cold bottom 
water south of the margin all indicate large heat transfers 
throughout the entire water column in the vicinity of the 
ice margin. The extreme temperature gradients associated 
with the presence of structure were found predominantly 
near the ice margin, but not far under the ice, indicative 
of strong downxvard mixing processes in the vicinity of the 
margin. Farther back under the ice, however, this powerful 
downward mixing of near-surface water appeared to be lacking. 
Therefore, ice keels probably are not the cause of mixing 
as was suggested by Corse (1974). Instead, the mixing must 
be due to a dynamic process which is peculiar to the ice 
margin. The process must be due to something distinctive 
of the ice margin, most likely the rapid melting of ice. 

This melting process and its relationship to the formation 
of mesostructure will be discussed in a subsequent section. 
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B. DENSITY CHARACTERISTICS 



The cross-sectional density distributions, plotted in 
sigma-t units, are shown in Figures 17 through 23. All 
crossings were typified by a high gradient layer in the 
upper 10 m, with localized pockets of low density water 
located at the surface which correspond to areas of melting 
ice. The density surfaces like the isotherms, displayed 
irregularity in the vicinity of the ice margin. However, 
the density surfaces, unlike the isotherms, displayed no 
inversions. The density at the surface was observed to be 
higher outside the margin and decreased near it. These 
phenomena were undoubtedly due to melting processes. The 
irregularity, high density gradients, and low density 
pockets did not appear farther back under the ice. The 
rapidly added melt water must have been mixed downward 
rather quickly in the vicinity of the margin so as to have 
decreased the high gradients that were present. The density' 
cross sections display the same general characteristics as 
do the temperature cross sections and there is evidence of 
dynamic activity in all cases. 

The deepest isopycnals, present well inside the ice 
margin, were destroyed by downward mixing of melt water 
near the margin. These isopycnals usually rose toward the 
interior of the ice field, forming a wedge-like feature. 

In Crossings 1, 2, 3, and 4 the wedge was pronounced. In 
Crossings 5 and 6 the wedge was less pronounced and did 
not exist in Crossing 7. This wedge was well correlated 
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Figure 17. Sigma-t Cross Section, Crossing 
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Figure 19. Sigma-t Cross Section, Crossing 
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Figure 20. Sigma- t Cross Section, Crossing 
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Figure 23. Sigma-t Cross Section, Crossing 



with the cold, bottom water found in the temperature 
data . 

The complex tilting of the density surfaces is evidence 
of lateral pressures that must be maintained by dynamic 
processes. Thus, stronger density structures should be 
indicative of strong dynamic processes and more severe 
mesostructure . Deep mesostructure was found to be more 
severe for those crossings with a well-defined wedge and 
thus gave support to 'this hypothesized relationship between 
deep mesostructure and strong dynamic processes. This 
relationship was further confirmed by an investigation of 
the dynamic heights associated with the marginal ice zone 
and is discussed in a subsequent section. 

Observations of the 26.0 and 26.5 sigma-t surfaces 
revealed that they had a pattern of variable spacing 
related to the ice edge. For Crossings 1, 2, 3, and 4 
(Figures 17-20) these isopycnals were more widely separated 
outside the ice margin, thus defining a zone of very low 
vertical density gradient. For Crossings 5, 6, and 7 
(Figures 21-23) this low density gradient zone was found 
well inside the ice edge, but in Crossing 5 it occurred 
between the 26.5 and 26.8 sigma-t surfaces. Inspection of 
the above figures indicates that most of the deep meso- 
structure is found within these zones of weak density 
gradient. This .is perhaps an expected phenomenon as 
relatively little energy is required to mix anomalous 
water downward in the water column when the density 
gradient is small. 
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It is important to remember at this point that a 
temperature or density cross section suggests that there 
is a transport of properties from left to right. This need 
not be the case; flow may be oblique or nearly parallel to 
the ice edge. This type of flow pattern may account for 
small degrees of tilt to the isopycnals within the ice 
or for weak dynamics and weak structure to exist. It may 
also cause a low rate of heat transport to the ice and 
consequently low melting rates. Considering these factors 
the heat balance across the ice margin was investigated. 

C. HEAT BALANCE ACROSS AN ICE MARGIN 

The disappearance of the warm surface water under the 
ice margin and the melting which occurred led to an analy- 
sis of the heat content of the water columns associated with 
the ice margin crossings. Specifically, it was desired to 
correlate the amount of heat dissipated at the margin with 
the melting of an equivalent amount of ice and to test the 
associated idea that water flows along the section toward 
the ice. If it did not, it might be an explanation for 
weak dynamic mixing processes and weak structure. 

The heat content of the water column at each station 
in the various crossings was calculated, using as a ref- 
erence water of 33.33°/oo and -1.8°C. This reference was 
selected since 33.33°/ 0 o is the salinity of freezing 
sea water at the lowest water temperature observed in the 
area. The resultant values were termed the "heat excess" 
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as defined by equation (1) : 



^D 



H = C 



y 



p(T-T R ) dz 



( 1 ) 



where 



_ ? 

H = Heat excess (cal-cm ) 



0^ = Specific heat (cal-gm 

_ 3 

p = Density (gm-cm ) 



-1 op-1- 



T = Temperature (°C) at depth z cm 
T r = Reference temperature (°C) 

D = Depth of the water column (cm) 
dz = Depth increment (cm) 

In order to examine the melting phenomenon in terms of 
the fresh water introduced from melting ice an "ice melt 
content" was also defined. It was defined as that weight 
of ice, with an initial salinity of 6.0°/ O o and a density 

_3 

of 0.94 gm cm , that must be transformed into melt water 
to produce the salinity structure present in the given 
column of water, referenced again to a water column having 
an initial salinity of 33.33°/oo- The ice melt content is 
defined by equation (2) : 



where : 



I 



r 



D 



P (■ 



R 






’R 



■q ) dz 
b I 



- 2 

I = Ice melt content (gm-cm ) 

S R = Reference salinity (ppt) 

S = Salinity (ppt) at depth z cm 
Sj = Salinity of ice (ppt) 



( 2 ) 
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Values of 0.956 cal gm ^ °C ^ for C and 100 cm for dz 

P 

were selected and the values for the remaining variables 
(T, D, p, and S) were obtained from the individual station 
data. Using these values the total integrated heat excess 
and ice melt content from the surface to the bottom were 
computed. The values from the surface to 35 m ivere selected 
for comparison purposes for the analysis. One may make 
nearly as good an argument for integrating to the bottom. 
However, the likelihood of incorporating an anomalous 
contribution from a variable depth limit led to the choice 
of a fixed limit. The results of these computations for 
all crossings are displayed in Table I. 

Ideally, with a flow of water directly across the ice 
margin, the heat excess and ice melt content data should 
present an inverse relationship. The heat excess should 
decrease and remain low after passing the margin as the 
ice melt increases, which in turn should remain high due to ' 
the melting process. Considering the heat excess and ice 
melt content data, many irregularities were present in each 
of the crossings. However, all crossings did display a 
general heat decrease and an increase in ice melt content 
close to the ice margin, except for Crossing 6 which 
displayed a decrease in ice melt content and heat excess. 
Crossing 2 is presented as an example (Figure 24) . It 
was close to the average of all the crossings, excluding 
Crossing 6, even to the extent of being less irregular than 
most. The data from the seven crossings revealed a decrease 
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Crossing Ice Melt Content (gm cm -2 ) * = Ice Margin 

Stations Heat Excess (cal cm- 2 ) 
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DISTANCE (NMI) 

Figure 24. Heat Excess and Ice Melt Content for Crossing 



_ 2 

in excess heat ranging from 4900-10300 cal cm across 

_ 2 

the ice margin. The average change was 6600 cal cm 

Equation (3a) describes the heat (AH) that is required 
to warm and melt ice of initial temperature T q and thick- 
ness h 



AH = [(T - T ) C + X T ]hp 

g o' p I J H 



(3a) 



where : 



Solving 



_ 2 

AH = Heat change (cam cm ) 

T q = Initial temperature (°C) 

T = Melting temperature (-.053S) (°C) 

S 

S = Salinity (ppt) 

Cp = Specific heat of fresh ice (cal gm“ lo C* 1 ) 

Xj = Latent heat of fusion of sea ice (cal gm ^) 
h = Ice thickness (cm) 

for the ice thickness, h, we arrive at equation (3b) 



- £ [(T 

P g 




c + 



p 




(3b) 



An ice cover of salinity 6.0 o / oo , initially at -2°C, 

- 7 . 

and with a density of 0.94 gm cm was selected as repre- 
sentative of the area of study. This selection resulted 

- 1 - 1 - 1 
in values of 0.5 cal gm °C for and 67 cal gm for 

Aj. Using the average change of heat content, 6600 cal cm 

across the margin and the above values, an ice thickness 

of 1.0 m would be melted. 
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The ice melt content, with one exception, displayed an 

. 2 

increase across the margins of 30-73 gm cm with an 

-2 

average value of 54 gm cm . The exception was Crossing 6 

which resulted in a negative value. The average value of 
_ 2 

54 gm cm ' when equated to a thickness of reference ice by 
equation (4) resulted in an equivalent ice thickness melted 
of 0.6 m: 

h = Alp" 1 (4) 

where: h = Thickness of ice melted (cm) 

_ 2 

AI = Change in ice excess (gm-cm ) 
p = Density of ice (gm-cm ) . 

The values of the thickness of ice melted computed 
from the heat excess and ice melt content agree favorably 
with the ice thickness of about 0.8 m predicted for the 
region of study by Wittman et al . (1963). This lends 
support to the assumptions that the dissipation of the heat 
across the ice margin can be accounted for by the melting 
of the ice cover and that the flow under the ice margin 
is generally northward (left to right in the sections) . 

The overall changes in heat excess and ice melt content 
across the margin were not apparently related to the severity 
of structure present (Table III). However, the rates at 
which they changed near the margin were possibly related 
to the severity of shallow mesostructure . Crossings 1 
and 2 had the least rapid heat and melt water changes in 
the vicinity of the margin and also displayed the least 
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severe shallow structure. Crossings 3, 4, and 5 had 
significantly more rapid heat exchanges at the margin and 
also had the more severe shallow structure. These factors 
suggest a correlation between the rate of heat exchange 
and the severity of shallow structure. 

In most cases the major portion of the heat exchange 
occurred at the ice margin; however, in Crossings 1 and 
2 it occurred outside the ice margin. This can be explained 
by a flow away from the ice subsequent to a period of 
melting or by a retreat of the ice front that was more 
rapid than the northward flow of water. 

In Crossings 1, 3, 4, 5, and 7 there was an irregular 
increase of ice melt content toward the north which is 
what one might expect if the melt water near the margin 
were carried under the ice by the current. Crossings 2 
and 6 displayed a gradual decrease of ice melt content 
as one progressed farther under the ice. Possible explana- ' 
tions for this difference are that the axis of flow of the 
northward current was not directly across the margin or 
an irregularity of flow occurred. Thus, it appears that 
the flow at the margins is 1 not as simple as was assumed, 
although it is generally northward. 

D. DYNAMIC HEIGHTS AND LOCAL CURRENTS 

The degree of disturbance in the density surface 
(Figures 17-23) and the presence of many factors which 
pointed to the existence of highly dynamic processes acting 
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in the immediate vicinity of the ice edge gave rise to an 
interest in the dynamic heights associated with melting at 
the marginal ice zone. Neshyba's model (1974) of the pro- 
grading ice cover during the freeze-up of the Bering Sea 
predicted the inverse phenomenon, a dynamic low under the 
ice margin. His model described the freezing process 
which rejected brine and mixed it downward through the entire 
water column resulting in higher density, a dynamic low, and 
a horizontal pressure gradient toward the ice under the 
freezing zone. Neshyba concluded that freezing rates were 
slow enough such that the pressure gradient must be in geo- 
strophic balance and must therefore correspond to a 
geostrophic flow to the east under the freezing zone. 1 1 is 
assumption that the oceanic system was otherwise static 
probably led to unrealistically high values of salinity 
and correspondingly high negative values of dynamic height. 

An extension of Neshyba's work to the melting season would 
suggest that the converse situation, a dynamic high, might 
be found under a melting ice zone and, with melting much 
more rapid than freezing, that the change in geopotential 
should be much greater. 

Dynamic heights for all crossings were computed utilizing 
the HYDRO program (W. R. Church Computer Center, NPS) and 
the results appear in Table II. Figure 25 presents a 
vertical cross-section of dynamic heights for Crossing 2 
which was selected to illustrate significant features. All 
crossings in general displayed local irregularities and 
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Table II. Dynamic Heights for Crossings 
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Figure 25. Dynamic Heights for Crossing 



all, with the exception of Crossing 6, displayed a hill 
(or half of a hill as in Crossings 1, 3, and 4) situated 
near the ice margin. Evidence of this dynamic hill was 
observed to a depth of about 10 m in Crossings 1, 2, 3, and 
4 and to about 15 m in Crossings 5 and 7. Thus the hill 
was generally noticeable in the near surface zone only. 

These results were somewhat as predicted. However, 
the hill was found to be a localized phenomenon extending 
over horizontal distances of up to 30 km across the margin 
and its effects were observed to only about 10 m of depth. 
This was different from Neshyba's freezing model which 
resulted in a general depression associated with the entire 
freezing zone which extended 100 km back under the ice, 
and whose effects were notable throughout the entire water 
column. The heights of these hills were of the order of 
2 dyn cm as compared to depressions of 2.7 dyn cm calcu- 
lated by Neshyba. The hills due to melting were expected 
to have been several times as large as the depressions due 
to freezing, but were not. This was probably because the 
water at the ice margin was not static as in Neshyba's 
model. It will be seen, however, that the magnitude of 
this hill is not the critical factor. 

The author's data revealed that the height of the hill 
under the ice apparently was not correlated with the severity 
of the mesostructure . However, the magnitude of the hori- 
zontal gradient of the dynamic height at the margin did 
appear to be related directly to the severity and presence 
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of deep mesostructure . Table III shows that Crossings 1, 

2, 3, and 7 had the highest dynamic height gradients and 
also the most severe deep structure, whereas Crossings 4 
and 5 had relatively small gradients and either lacked 
or had very little deep structure. 

The gradient observed in the dynamic heights must 
result in a lateral pressure gradient. This gradient is 
probably of recent origin due to the rapid retreat of the 
ice margin; thus geostrophic balance is probably not achieved 
and the pressure gradient is most likely balanced against 
inertial acceleration or internal friction. In either case 
the pressure gradient leads to a velocity shear in the 
vertical plane normal to the ice margin. Because the 
pressure gradients on the two opposing sides of the dynamic 
hill are oppositely directed, there will be a diminution 
of flow near the surface on the seaward side of the dynamic 
hill and an acceleration of flow on the other. Continuity 
of flow then would demand that the shear at the surface 
be compensated by a reverse shear at deeper depths under 
the dynamic high. It is believed that this downward directed 
flow may be the agency which caused the mixing at depth 
which was observed near the ice margin. 

Both of the possible opposing forces of friction and 
inertia were found to be of the same order of magnitude, 
depending upon one's choice of viscosity and characteristic 
velocity, but viscosity is likely to exert the largest 
effect. Therefore, to explore the magnitudes of the 



57 



4 -> 

o a cd 

*H CD rJ^ 

& *H PJ 

aj 

C oj P aj 
XH aj^ 
Q C3 



a > 

H £ 

bOP4 
P 



CO 

to 

o a 

P *H /- 
O bOrM 
aj P I 
Pi aj g 
•HPS 

0 CD CD 
+J U W 
PI l-H 

_ o < 

Pi ej CD o 
O iH 

a p> +j 

HH i— I 
CD 



CD 

a 



to 

aj 

0 



CM 

0 PI I 

pp *h g 
P P W)U 
p 

0 0 to aj p- 1 
P E to S oj 
o u 
p 0 
O O <*> 
aj h o 



0 
to 
aj aj 



O 

0 a 

Q *H 



0 xt 
p 0 

0 P5 r-H a 
tn P H -H 
O aj 0 m 
^0^0 
P-* Q 



0 

0 X P 
P P 0 
3 > 

P to lO 0 
U Pi f CO 
p 0O t 
P P> ' 0 



+-> a 

CO IH 



PI 

o 



bi) 

PJ 

•H 

to 

to 

o 

p 

u 



p. 

0 

0 

o 



aj 

X 

CO 



CXI 


Cb 


LO 


o 


vO 


oo 


r^ 


CXI 


to 


CXI 


1— 1 


rH 


o 


cxi 


o 


o 


o 


o 


O 


iH 


o 


o 


o 


o 


o 


o 


O 


o 












1 




o 


to 


vO 


o 


H* 


H* 


oo 


to 


t'- 


vO 


LO 


vO 


O 


to 












tH 

l 




r^ 


vO 


H* 


CXI 


cn 


cr> 


to 


• 


• 


• 


• 


• 


• 


• 


LO 


vD 


vO 


vO 


H- 


LO 


o 














iH 








to 


to 


to 


to 


o 


O 


O 


0 


0 


0 


0 




'Z 




>-< 


>h 


>H 


>< 


LO 


LO 


LO 


rH 


rH 


o 


(XJ 


iH 


rH 


LO 


CXI 


to 


to 


OO 


pH 


CXI 


to 


H* 


LO 


vO 

1 


r-- 



58 



Table III. Comparison of Crossings. 



effects, approximate calculations were made based on the 
assumption that the lateral pressure gradient was balanced 
by the frictional shear of eddy viscosity. A model was 
developed assuming a current, u, outside the ice, flowing 
in the plane of the pressure gradient and uniform in 
velocity throughout the upper warm layer. 

The step-like profile of Figure 26 (solid line) was 
designed for the current velocity structure hypothesized 
outside the ice margin based on observed currents from 
other investigations. Equation (5) was used to describe 
the balance of frictional shear with the horizontal pressure 
gradient : 

A a 2 u = AD , . 

p 9z 2 - AL ^ 

A = Coefficient of eddy voscosity (gm-cm ^ sec ^) 

- 3 

p = Density (gm-cm ) 
u = Current speed (cm-sec ^) 
z = Depth (cm) 



AD . 
AT " 



Horizontal 



gradient of dynamic height 



(dyn m-km 1 ) 



Equation (5) was solved for u with the conditions that the 

dynamic effect was limited to 10 m and that z = 10 m, 

= 0 and u = 15 cm sec Specification of this velocity 

was not necessary to the solution, but it did influence to 

3 2 U 

some extent the choice of A. ^ — ? was assumed to be constant 

a Z 

which required that p and A also be constants. 
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Figure 26. Effect of the Dynamic Hill 
on the Current Profile. 
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The values of the horizontal gradient of dynamic 
height, which were observed for the crossings, ranged 
between 0.0010 and 0.0039 dyn m km ^ with an average of 
0.0026 dyn m km ^ for all crossings except Crossing 6 
(Table III). Values for density and depth were selected as 

_3 

1.02 gm cm and 10 m, respectively. Selecting a value for 

the coefficient of eddy viscosity (A) was more difficult. 

Eddy viscosity coefficients have an extreme range of values 

which vary with the scale and rate of the phenomenon in 

question. Neuman (1966) shows values of from 1 to 3000 
- 1 - 1 

gm cm sec . Campbell (1965), however, arrived at a 
value of 55 gm cm ^ sec ^ for a numerical model of the ice 
packs of the Arctic Ocean, in which a relative velocity of 
13 cm sec ^ existed between an ice floe and the sea water. 
Figure 27 shows how a variation in A affects the magnitude 
of the total velocity diminution experienced at the surface 
due to the shear forces opposed to the flow. Varying A 
from 50 to 500 gm cm ^ sec ^ reduces the current from 23.5 
to 2.35 cm sec ^ which is of the proper order of magnitude 
for the currents assumed for the region. An intermediate 
value of 100 gm cm sec resulted in a decrease of the 
velocity at the surface of 13 cm sec ^ , a fortuitous agree- 
ment with Campbell. 

The application of this velocity shear in the upper 10 m 
of the water column and the enforcement of continuity of 
flow upon the initial current profile of Figure 26 resulted 
in the second profile (dashed line) of Figure 26, the 
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VELOCITY (CM sec 1 ) 




Figure 27. Effect of Eddy Viscosity on 
Current Velocity Diminution. 
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profile for the flow beneath the ice margin. This profile 
displays a descent of the core of the current as it 
encounters the margin and its effects. The compensating 
flow at depth is pictorial; it did not seem feasible to 
calculate it in a more meaningful way. This restriction to 
the flow, the resultant downward motion, and the accelera- 
tion of flow at greater depths could provide a mechanism 
for the mixing process which occurred. 

It is evident that where dynamic height gradients are 
large, the vertical shear must be large and require greater 
and probably deeper compensating flows at depth. As pre- 
viously discussed, these regions will most likely be areas 
of significant deep structure. 
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IV. CONCLUSIONS 



A. SUMMARY OF RESULTS 

The distribution of temperature, density, heat, and 
dynamic heights confirmed that the phenomenon of temperature 
mesos tructure near Arctic ice margins was related to the 
presence of an ice margin and the associated melting pro- 
cesses. Mesos tructure was characterized by large tempera- 
ture gradients occurring simultaneously with low density 
gradients. It was present in a relatively broad zone, 

18 to 65 km wide, across the margins. It varied consider- 
ably over short distances, with individual elements 
generally being traceable for only a few kilometers. 

The assumption of a generally northward flow for the 
area of study was supported by the heat excess and ice 
melt content calculations. A decrease in heat content and 
a corresponding increase in the ice melt content of the 
water column was observed as one crossed the ice margin. 
These changes agreed reasonably well enough with the heat 
required and the melt water that would have been formed by 
the melting of the existing ice cover to support the 
assumption of a generally northward flow. The flow, 
however, displayed the possibility of not always flowing 
directly across the margins. Crossings 2 and 6 displayed 
both a heat and an ice melt content decrease as one 
progressed farther under the ice. This could have 



64 



indicated a flow approaching the margin obliquely and thus 
carrying some of the melt water away. The many irregulari- 
ties also present in both the heat and ice melt results 
suggest the presence of complex flows near the margins. 

Large heat exchanges and downward mixing throughout 
the entire water column were evident near the ice margins. 

The warm surface flow was cooled at the surface by melting 
to form the nose feature. The cold, dense bottom water, 
not present south of the margin, was dissipated by downward 
mixing in the water column creating a high density wedge- 
like feature. The broad zones of low vertical density 
gradient resulting from this mixing were found only near 
the margin and not far under the ice. Thus, the creation 
of mesostructure by the turbulent mixing by ice keels 
as suggested by Corse (1974) is probably a doubtful 
mechanism due to the lack of evidence of mixing farther 
under the ice. 

The area near the ice margin was characterized by large 
temperature gradients along isopycnals and by the irregularity 
and tilting of the isopycnals. These irregularities and 
tilting must result in lateral pressure gradients that 
evidently are balanced by complex dynamic flows in the 
marginal area. 

Further investigation led to the discovery of a dynamic 
hill which extended about 30 km across the margin. The hill 
was due to the melting process and was notable to a depth 
of about 10 m. It was interesting that the dynamic hills 
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were about the same height as the depression computed by 
Neshyba (1974) for the freezing process. If his numbers 
are correct, the dynamic hill for the melting case should 
be many times larger since melting occurred much more 
rapidly than freezing and was concentrated over a much 
narrower width. However, it appears that fresh water was 
lost from the area by circulation, thus greatly diminishing 
the freshening effect. 

The magnitude of the dynamic hill was found to be of 
lesser importance than the gradient of the hill near the 
ice margin. The lateral pressure gradient which resulted 
from the horizontal gradient in dynamic height was assumed 
to be balanced by internal frictional stresses, and when 
approximate calculations of the effect of this balance were 
made, a velocity shear at the surface resulted. This velocity 
shear was applied to a simple velocity profile which resulted 
in a deepening of the flow axis and an acceleration of the 
flow at depth. This result provides a plausible explanation 
for the source of the downward mixing that is observed 
at the margins. This hypothesis is also supported by the 
correlation of larger horizontal dynamic gradients at the 
surface with a greater sloping of the deep isopycnals and 
the presence of more severe deep mesos tructure . The 
greater dynamic forces at the surface should result in 
larger velocity shears and thus stronger effects at depth. 

Shallow and deep mesos tructure are believed to be 
related to different degrees of mixing action. Shallow 
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structure was correlated with the presence of the nose 
feature while deep structure was related to large dynamic 
slopes. Deep structure was usually within the density zone 
of 26.0-26.5 sigma-t units which was found below 20 m 
depth. The deep structure was located in the vertically 
broadened areas of this density zone. The low density 
gradient of this broadened area implies mixing processes, 
which one might expect from the large gradients in dynamic 
height. This area of low vertical density gradient also 
was one where interleaving might easily occur due to the 
small amounts of energy required to lift a parcel of 
water. Deep structure appeared to be more severe when the 
nose feature was not well defined. Thus it would appear 
that the higher dynamic slopes at the surface cause more 
severe velocity shears. The resulting stronger mixing 
effects can destroy the nose feature and thus explain the 
usual absence of the nose when deep structure is strong. 

Shallow mesostructure was related to the nose feature 
and was generally found above the 26.0 sigma-t contour and 
thus usually shallower than 20 m depth. It was observed 
at the boundaries of the nose feature where water types of 
similar density and different temperatures were available 
to interleave. Shallow structure, with two exceptions, 
was found to be more severe in the cases in which deep 
structure was weak or absent. 
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B. RECOMMENDATIONS FOR FUTURE INVESTIGATIONS 

As a result of this study it is evident that additional 
data are necessary in future investigations. The probable 
variations in local flows near the ice margin make it 
apparent that current vector data are needed for the stations 
taken near the margins. Data also need to be collected far- 
ther into the ice and more exact data on ice concentration 
and thickness would be helpful. Time series data on the STD 
parameters and currents, taken simultaneously, would also 
be useful in determining the time variance of the phenomena. 
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APPENDIX A 



STD DATA FOR STATIONS ANALYZED 

The data for each station of the ice crossings analyzed 
follows. Four stations appear on each page and the respec- 
tive station number appears in the lower part of each plot. 
The scales for the horizontal axes are located at the top 
of each column and the symbols are as follows: 

ST, sigma-t 

SV, sound velocity (m sec 

S, salinity (°/ 0 o) 

T, temperature (°C) 

The symbol for each of the above parameters is at the bottom 
of each trace and an S, T, or V (salinity, temperature, or 
sound velocity) appears at the top to help distinguish indi- 
vidual traces. The gap usually located at about 10 depth 

i 

on each trace is due to the requirement for two lowerings 
of the STD. 
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